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Abstract: We present an extensive study of non-minimal flavour violation in the squark 
sector in the framework of the Minimal Supersymmetric Standard Model. We investigate 
the effects of multiple non-vanishing flavour-violating elements in the squark mass matrices 
by means of a Markov Chain Monte Carlo scanning technique and identify parameter 
combinations that are favoured by both current data and theoretical constraints. We 
then detail the resulting distributions of the flavour-conserving and flavour-violating model 
parameters. Based on this analysis, we propose a set of benchmark scenarios relevant for 
future studies of non-minimal flavour violation in the Minimal Supersymmetric Standard 
Model. 
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1 Introduction 

Among all possible extensions to the Standard Model of particle physics, weak scale super- 
symmetry [1, 2] remains one of the most popular and best studied options. The quest for 
the superpartners of the Standard Model degrees of freedom is one of the hot topics of the 
current high-energy physics experimental programme and many search channels are hence 
investigated at the Large Hadron Collider (LHC) [3, 4]. Since no signal of supersymme¬ 
try has been found so far, the results have been interpreted either in terms of limits on 
specific setups like constrained versions of the Minimal Supersymmetric Standard Model 
(MSSM) or in terms of simplified model spectra inspired by the MSSM. As a result, ei¬ 
ther supersymmetric particles are constrained to reside at scales that are not reachable 
in proton-proton collisions at a centre-of-mass energy of 8 TeV, or the spectrum must 
present specific properties that allow the superpartners to evade detection as in the case 
of compressed supersymmetric spectra or non-minimal realizations of supersymmetry. In 
this work, we follow this latter guiding principle and explore to which extent deviations 
from the minimal flavour-violation paradigm [5-7] are allowed by current data. 

In minimally flavour-violating supersymmetry, all the flavour properties of the model 
stem from the diagonalization of the Yukawa matrices that yields different gauge and mass 
bases for the (s)quarks and (s)leptons. Flavour violation is thus entirely encompassed 
in the CKM and PMNS matrices. There is however no theoretical motivation for the 
flavour structure of a supersymmetric model to be the same as in the Standard Model. For 
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instance, when supersymmetry is embedded in a Grand Unified framework, new sources of 
flavour violation could be allowed [8] . The soft supersymmetry-breaking mass and trilinear 
coupling matrices of the sfermions could therefore comprise non-diagonal flavour-violating 
entries that are not related to the CKM and PMNS matrices. We adopt such a setup, which 
is referred to as non-minimally flavour-violating (NMFV) supersymmetry, and study the 
impact of these additional flavour-violating soft terms. 

In recent years, the consequences of non-minimal flavour violation in the squark sector 
have been investigated in various areas. NMFV effects on low-energy observables such as 
rare decays (see, e.g. Ref. [9] and references therein) or oblique parameters [10] have been 
considered, and the potential signatures at the LHC have been investigated [11-24]. More 
recently, the existing constraints on possible non-vanishing flavour-mixing parameters have 
been updated [25, 26]. These results have been derived under the restriction that only few 
off-diagonal elements of the squark mass matrices are non-zero, and that at most two of 
them are varied at the same time. One would however generally expect that several of 
the flavour-violating entries could be non-vanishing, especially if the flavour structure is 
generated by some mechanism at a higher scale. Consequently, a comprehensive study of 
the most general NMFV configuration of the MSSM, where all flavour-violating Lagrangian 
terms are taken into account and confronted to current data and theoretical constraints, 
is in order. A first step in this direction is achieved with this work. 

We consider the most general mixings between second and third generation squarks. 
Any non-CKM induced mixing with the first generation is ignored as a result of constraints 
imposed by kaon data [27]. Choosing a phenomenological approach, we model the flavour- 
violating effects under investigation by a set of 19 free parameters defined at the TeV scale 
and identify, by means of a Markov Chain Monte Carlo parameter scanning technique, 
the regions of the parameter space that are favoured in light of current data. The rest of 
this paper is organised as follows. We first review in Section 2 the squark sector of the 
MSSM when NMFV is allowed and present our parameterisation of the effects under study. 
In Section 3, we describe the machinery that is employed to explore the 19-dimensional 
NMFV MSSM parameter space and present the experimental constraints that are imposed. 
Our results are discussed in Section 4, we propose NMFV MSSM benchmark scenarios for 
the LHC Run H in Section 5 and our conclusions are given in Section 6. 


2 The squark sector with general flavour mixing 


Starting from the most general MSSM Lagrangian, the super-CKM basis is defined by 
rotating the quark and squark superfields in flavour space in a way in which the quark mass 
matrices and nid are diagonal. Squark and quark flavours are thus aligned, although 
the squark mass matrices are not necessarily diagonal. In the {uL,CL,iL,uji,CR,iji) and 
{cIl, SL,bL, dR, SR,bR) bases, the up-type and down-type squark mass matrices A4| and 
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are given by 
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respectively. In these notations, we have introduced the soft supersymmetry-breaking 
squark mass matrices and M"^ for left-handed, up-type right-handed and down- 

type right-handed squarks respectively, as well as the matrices and that embed the 
trilinear soft interactions of the up-type and down-type squarks with the Higgs sector. 
While these five matrices are defined to be flavour diagonal in usual constrained versions of 
the MSSM, our NMFV framework allows them to be general and possibly flavour-violating. 
Moreover, Vckm stands for the CKM matrix, /r denotes the superpotential Higgs(ino) mass 
parameter and tan/3 = (((( is the ratio of the vacuum expectation values of the neutral 
components of the two Higgs doublets. Finally, the squark mass matrices also include 
(flavour-diagonal) H-term contributions 


Dq,L = ■'TT'Z {iq — Cg siu^ 9w) COS 2/3 and Dq^R = m\eq sin^ 9w cos 2/3 , (2.2) 


where mz is the Z-boson mass, 9w is the weak mixing angle and Cg and Ig (with q = u,d) 
are the electric charge and the weak isospin quantum numbers of the (s)quarks. 

In order to reduce the number of supersymmetric input parameters, we assume that 
the first two generations of squarks are degenerate so that the (flavour-conserving) soft 


masses are determined by six free parameters, 

(^q)ii = (^< 3)22 = ^Qi,2 ’ 

{^0)33 = ’ 
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to 

{^0)33 ^ ’ 

(2.3) 

(-^£>)ll “ {^0)22 — -^£>1,2 ’ 

{^0)33 ^ ■ 



Moreover, we define the diagonal components of the trilinear couplings Tg relatively to the 
Yukawa matrices Yq, 

{Tq)ii = {Aq)^^{Yq)ii . (2.4) 

We then neglect the first and second generation Yukawa couplings so that only the trilinear 
coupling parameters related to third generation squarks are considered as free parameters. 
We take them equal for simplicity, so that we have 


(^n )33 = At , {Ad )33 — A}, and At — A}, — Af . (2.5) 

We now turn to the off-diagonal elements of the squark mass matrices. In order to 
be compliant with kaon data, we ignore any mixing involving one of the first generation 
squarks [27]. Next, following standard prescriptions [28], we normalise the remaining non¬ 
diagonal entries of the squared squark mass matrices with respect to the diagonal ones and 
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parameterise all considered NMFV effects by seven dimensionless quantities 
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The physical squark states Ui and di (with f = 1,..., 6) are obtained by diagonalizing 
the squared squark mass matrices M.\ and A^| according to 

diag(m?^, m|,..., m|) = nqMln\ for q = u,d . (2.7) 

By convention the mass eigenstates are taken ordered such that < • • • < . The 

6x6 rotation matrices TZa and carry the information about the flavour decomposition 
of the squarks, 


— TZu (ul Cl ii Ur cr Ir^ , 

(^di d2 ds di ds de) = (di sr Ir (Ir sr 

and their different entries directly appear in couplings of the squarks to the other particles 
(see, e.g. Refs. [11, 18]). 

In addition, the gaugino sector is chosen to be determined by a single parameter, the 
bino mass Mi. The wino and gluino tree-level masses M 2 and M 3 are then obtained by 
making use of a relation inspired by Grand-Unified theories. 

Ml = -M 2 = -M 3 . (2.9) 

2 6 

The slepton sector is defined in a flavour-conserving fashion, so that the soft terms consist 
of three (diagonal) mass parameters that we set to a common value 

(M|^)ll = (M ^)22 = (M ^)33 = M^ , ( 2 . 10 ) 

and the slepton trilinear coupling matrix to the Higgs sector Ti contains a single non-zero 
entry, 

(r^)33 = TrH, = YrAf . (2.11) 

All flavour-conserving trilinear sfermion interactions with the Higgs bosons are conse¬ 
quently driven by a single input parameter Aj. The model description is completed by the 
definition of the Higgs sector that is parameterised in terms of the // parameter, tan /3 and 
the pole mass of the pseudoscalar Higgs boson niA- 
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Parameter 

Value 

Parameter 

Value 

a~^{mz) 

127.934 

rrie 

510.9989 keV 

mz 

91.1876 GeV 

rrif. 

105.6583 MeV 

Gf 

1.16637 X 10-5 GeV-2 

mr 

1.77699 GeV 

m^®(mc) 

1.25 GeV 

^CKM 

0.2272 

mf^(2 GeV) 

120 MeV 

^CKM 

0.818 

mp(2 GeV) 

3 MeV 

^CKM 

0.221 

mfS (2 GeV) 

7 MeV 

^CKM 

0.34 


Table 1. Standard Model sector of our NMFV MSSM parameter space. 

3 Setup and constraints 
3.1 Numerical setup 

In the previous section, we have defined a simplified parameterisation of general NMFV 
MSSM scenarios in terms of 16 soft supersymmetry-breaking parameters to which we have 
supplemented three parameters related to the Higgs sector. Turning to the Standard Model 
sector, the QCD interaction strength is computed from the value of the strong coupling 
constant at the Z-pole as{mz), while we choose as three independent electroweak inputs 
the electromagnetic coupling constant evaluated at the Z-pole a{mz), the Fermi constant 
Gf and the Z-boson mass mz- The fermion sector is defined by the pole mass of the top 
quark the MS mass of the bottom (charm) quark rub (me) evaluated at the mb {me) 

scale and the MS masses of the three lightest quarks evaluated at a scale of 2 GeV. Finally, 
we include in our parameterisation the masses of the electron {me), the muon (m^) and 

the tau {mr) and we calculate the CKM matrix using the Wolfenstein parameters 
TCKM^ pCKM ^CKM_ 

For our exploration of the NMFV MSSM, we start by fixing the Standard Model pa¬ 
rameters to the values provided in the review of the Particle Data Group [29], as shown in 
Table 1. We then allow as{mz), m^°*® and m^^{mb) to vary according to Gaussian profiles 
with their measured central values and errors taken as means and widths, as shown in Ta¬ 
ble 2. We vary randomly all the 19 supersymmetric parameters in the intervals presented 
in this table and finally calculate the resulting mass spectrum of the supersymmetric par¬ 
ticles at the one-loop level using the publicly available spectrum generator SPheno [30, 31]. 
We can in this way widely cover the regions of the parameter space where the electroweak 
symmetry is successfully broken (with, e.g. sub-TeV values for the /r parameter) and whose 
signatures are in principle observable at the LHC within the next few years (i.e. with not 
dramatically large soft sfermion masses). Moreover, the trilinear coupling parameter Af 
and the d-parameters have been chosen to prevent all off-diagonal elements of the squark 
mass matrices from being too large, so that tachyonic mass-eigenstates are avoided. 

In order to explore the 22-dimensional parameter space summarized in Table 2, we 
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Parameter 

Scanned range 

Parameter 

Scanned range 

as{mz) 

AA(0.1184, 0.0007) 

tan/3 

[10, 50] 

pole 

AA(173.3,1.3928) GeV 


[100, 850] GeV 

mb{mb) 

AA(4.19,0.12) GeV 

rriA 

[100, 1600] GeV 

Mqu. 

[300, 3500] GeV 

Ml 

[100, 1600] GeV 

Mqs 

[100, 3500] GeV 

Mi 

[100, 3500] GeV 


[300, 3500] GeV 

^LL 

[-0.8, 0.8] 

Mus 

[100, 3500] GeV 

^RR 

[-0.8, 0.8] 


[300, 3500] GeV 

°RR 

[-0.8, 0.8] 


[100, 3500] GeV 

°LR 

[-0.5, 0.5] 

A 

[-10000, 10000] GeV 

su 

°RL 

[-0.5, 0.5] 


or \Af \ < 4max{Mq,M^} 

^LR 

[-0.05, 0.05] 



^RL 

[-0.05, 0.05] 


Table 2. Supersymmetric and Higgs sectors of our NMFV MSSM parameter space, as well as 
varying Standard Model parameters. denotes a Gaussian profile of mean ^ and width a. 


rely on a Markov Chain Monte Carlo (MCMC) scanning technique [32-34] and impose on 
each of the studied setups a set of constraints that is described in the next subsection. 
In this scanning procedure, a given point is accepted or rejected based on the comparison 
of the products of likelihoods of this point with that of the previous point, where each of 
the likelihoods is associated with a specific constraint accounting for measurements and 
theoretical predictions in the NMFV MSSM framework. 

3.2 Indirect constraints on general squark mixing 

The masses and flavour-violating mixings of the superpartners can be indirectly probed by 
numerous flavour physics constraints, the anomalous moment of the muon as well as by 
the properties of the recently discovered Standard-Model-like Higgs boson. In our MCMC 
scanning procedure, we additionally impose the lightest superpartner to be the lightest 
neutralino, so that it could be a phenomenologically viable dark matter candidate. We 
dedicate the rest of this section to a brief description of all observables that have been 
considered in the scan and that are summarised in Table 3. 

Non-minimal flavour-violating squark mixing involving third generation squarks is by 
construction very sensitive to constraints arising from H-physics observables. In particular, 
H-meson rare decays and oscillations are expected to play an important role as the Standard 
Model contributions are loop-suppressed. Although we only consider squark mixing be¬ 
tween the second and third generations, we also include constraints arising from observables 
related to the kaon sector. Even if not present at the scale at which we calculate the super- 
symmetric spectrum (i.e. the electroweak symmetry breaking scale), squark mixings with 
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Observable 

Experimental result 

Likelihood function 

BR{B ^ V,7) 

BR{Bs —i /i/i) 

BR{B GeV^ 

AFB(H —)• 1 6] GeV^ 

BR{B -j- V5/X^)g2g[i^6] GeV^ 
BR(H —)■ Vs/i^)g2>;^4 4 QgV2 
BR(H„ tv)/BR{Bu rz2)sM 
AMb^ 

(3.43 ± 0.22) X lO-'^ [35] 

(2.8 ±0.7) X 10“® [36] 

(1.7 ±0.31) X 10"'^ [37] 
(-0.075 ± 0.036) X 10"^ [38] 
(0.66 ± 0.88) X 10-® [39] 

(0.60 ± 0.31) X 10“® [39] 

1.04 ±0.34 [29] 

(17.719 ±3.300) ps"^ [29] 

Gaussian 

Gaussian 

Gaussian 

Gaussian 

Gaussian 

Gaussian 

Gaussian 

Gaussian 

BR{K^ TT^vv) 

BR(Ar+ ^ TT+vv) 

(2.228 ± 0.29**") X 10"^ [29] 

< 2.6 X 10“® [29] 

1.73t|;J^ X 10“^° [29] 

Gaussian 

1 if yes, 0 if no 

Two-sided Gaussian 

Aa^ 

(26.1 ± 12.8) X 10-*® [e+e-] [29] 

Gaussian 

rrih 

125.5 ± 2.5 GeV [40, 41] 

1 if yes, 0 if no 

Lightest supersymmetric particle 

Lightest neutralino 

1 if yes, 0 if no 


Table 3. Experimental constraints imposed in our scan of the NMFV MSSM parameter space. 


the first generation are induced by the non-vanishing CKM matrix and renormalisation- 
group running so that kaon physics observables (calculated at a different scale) are also 
relevant for extracting constraints on the NMFV MSSM parameter space. 

We focus on the branching ratios associated with the rare B —)• Xgj, B —?■ K* 

B — )• XsHfj, and Bu —t tv decays, as well as on the forward-backward asymmetry (AFB) 
arising in —?• decays. The associated predictions are calculated with the Superlso 

package [42, 43]. In addition, we compute the neutral B-meson mass difference AMb^, 
the branching ratio associated with the Bg —)• —)• t^vv and K~^ —)• t^vv 

decays and the kaon parameter sk with the SPheno code [30, 31]. We furthermore employ 
SPheno for the estimation of the supersymmetric contributions to the anomalous magnetic 
moment of the muon and for a calculation of the lightest Higgs boson mass mh. The 
data transfer between SPheno and Superlso is achieved through the Flavour Les Houches 
Accord standard [44], and the Wilson coefficients for all hadronic observables are calculated 
in SPheno (at the scale Q = 160 GeV) and Superlso (at a scale Q = mw) from the values 
of the running coupling constants and the supersymmetric masses and parameters that 
have been evaluated with SPheno. 

We now briefly collect all references where the formulas that have been employed for 
the calculation of the considered observables can be found, and we indicate which NMFV 

^Imposing predictions for to agree with the related measured values leads to a preference for a lighter 
slepton mass spectrum. This can indirectly imply constraints on the NMFV MSSM parameters via the 
flavour observables that involve sleptons. 
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MSSM parameters are mainly constrained by each of these observables. The calculation of 
the B —7- Xg'y branching ratio is mainly based on the results of Refs. [45-47], while those 
of the B —7- and B —)• K*fx~ branching ratios respectively follow Refs. [48-51] 

and Refs. [52-57]. All three decays are sensitive to the left-left and to the left-right squark 
mixing parameters. 

In the case of the Bg —)• /ibranching ratio, the formulas of Ref. [58] have been 
used. For large values of tan /3, the pseudoscalar Higgs boson contribution gives a size¬ 
able deviation from the Standard Model expectation so that when non-minimal flavour 
violation in the squark sector is allowed, this observable mainly restricts left-right mix¬ 
ing parameters [59]. Additionally, it is also sensitive to Sll and when the gluino is 
not too heavy [60]. For the R-meson oscillation parameter AMb^, we use the formulas of 
Refs. [59, 61] with the hadronic parameters = —0.71, = —0.9, = —0.37 and 

pSLL _ _Q 72 . The NMFV contributions are mainly sensitive to the ^lr^rl 

SpR^BL products, the relative suppression and enhancement of their contributions being 
driven by the ratio of the chargino over the gluino mass [14, 59]. 

In the kaon sector, the observable is estimated by combining the formulas of 
Refs. [59, 62], the loop-contributions being evaluated with rju = 0.5, r]ct = 0.47 and 
r]cc = 1.44 [63]. In addition, we fix all hadronic parameters at the scale Q = 2 GeV 
as BY^^ = 0.61, = 0.76, BI^^ = 0.51, = 0.96 and B^^ = 1.2 [62], and 

we set the decay constant Jk to 155.8 MeV. The quantity ex is not directly sensitive 
to a single NMFV MSSM parameter but will allow us to constrain Sll, 13 ^bb 23 
5ll,23^bb13 products (recalling that first generation squark mixings are generated by 
renormalisation-group running). On different lines, the branching ratios associated with 
the rare K~^ —)• and Kb —t- decays are calculated from the formulas given in 

Ref. [64] with bl = 2.1310“^^, = 5.1610“^^ and Pc = 0.39. These observables mainly 

constrain the product 13 ^lr 23 higher-order combinations of (5-parameters 

that in particular appear in gluino/down-type squark box-contributions [64, 65]. 

In all the calculations described above, we have used the results of Ref. [66] for cal¬ 
culating chirally-enhanced interaction strengths that include, e.g. the resummation of 
loop-induced holomorphic coupling effects when tan/3 and/or the sfermion-Higgs trilinear 
couplings are large. 

As we allow for relatively light sleptons, charginos and neutralinos, we calculate the 
supersymmetric contributions to the anomalous magnetic moment of the muon by using 
the formulas of Ref. [67]. The related impact in terms of constraints on our NMFV MSSM 
parameter space depends on the higgsino/gaugino nature of the lighter charginos and neu¬ 
tralinos. 

Finally, the calculation of the Higgs boson mass includes the complete one-loop con¬ 
tribution that embeds all possible flavour structures and that is obtained by extending 
the formulas of Ref. [68]. For the two-loop corrections, we have made use of the formulas 
of Refs. [69-74] where generation mixing is neglected so that only third generation mass 
parameters in the super-CKM basis are used as input parameters. Although flavour ef¬ 
fects can shift the Higgs mass by a few GeV at the one-loop level, in particular when the 
product SpB^RL large [21, 26, 75], the two-loop effects are expected to be of one order 
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of magnitude smaller so that ignoring the associated flavour mixing is expected to be a 
reasonable approximation. 

4 Results and discussion 

The analysis of the results of the Markov Chain Monte Carlo scan presented in Section 2 
gives us information on the regions of parameter space that are favoured by the experi¬ 
mental data shown in Table 3. The influence of a specific experimental result on a given 
parameter can be studied by comparing its theoretical prior distribution to the posterior 
one that is derived after imposing the related constraint. The prior distributions of all 
parameters are obtained from a uniform random scan in which we ignore scenarios that 
exhibit tachyons, where the electroweak symmetry is not successfully broken and where 
the lightest neutralino is not the lightest supersymmetric particle. We hence include about 
1.5 X 10® theoretically accepted setups. The posterior distributions are then computed on 
the basis of our MCMC scan in which all the experimental constraints, except the one on 
the Higgs boson mass, are imposed. This scan consists of 100 chains of 6000 scenarios in 
which the first 900 ones (the burn-in length) are removed. The constraint on the Higgs bo¬ 
son mass is eventually imposed and the hnal posterior distributions include about 100 000 
points. To estimate the importance of each observable separately, we have run a separate 
MCMC scan consisting of 100 chains of 2000 scenarios for each observable. After removing 
the burn-in length, 170 000 points remain. For each scan, the convergence test of Gelman 
and Rubin has been verified [76]. 

Although we are mainly interested in the non-minimally flavour-violating parameters 
defined in Eq. (2.6), we first also discuss for completeness the flavour-conserving parameters 
of our model description. 

4.1 Flavour-conserving parameters 

We start the discussion with the twelve flavour-conserving parameters of our NMFV MSSM 
description. Figure 1 shows their probability density distributions over the respective 
parameter ranges. In each panel, we display the theoretical prior (yellow area) as well as 
the posterior distribution (solid line), which shows the impact of all constraints given in 
Table 3 together. 

The prior distribution of the gaugino mass parameter Mi is centred at relatively low 
values of Mi ~ 400 GeV and may reach values ranging up to about 1000 GeV. When im¬ 
posing all considered experimental constraints, the distribution is shifted by about 100 GeV 
to higher values. This feature can be traced to the chargino contributions to the Bg —)• P/U 
branching ratio and to the neutral H-meson mass difference AM^^, as we have fixed the 
ratios of the gaugino mass parameters Mi and M 2 so that chargino effects are connected 
to Ml. 

For the trilinear coupling parameter Aj, the prior distribution is centred around 
zero. Large values of Af are indeed often rejected since they can induce a large left- 
right squark mixing implying tachyonic states. Imposing the experimental constraints 
drastically changes the shape of the distribution and leads to two peaks corresponding 
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Figure 1. The one-dimensional prior (yellow histogram) and posterior (violet curve) distributions 
of the parameters of our NMFV MSSM description. The prior only incorporates theoretical inputs 
while the posterior distribution shows the impact of all experimental observations listed in Table 3. 

















to \Af\ r\j 3000 GeV. This feature is induced by the Higgs boson mass requirement that 
necessitates, in order to be satisfied, a relatively large splitting of the masses of the squarks 
exhibiting the largest stop component rriq-^ and rriq^. More precisely, the flavour-conserving 
formula for the leading contributions to the Higgs mass. 


2 2 2 na < 

cos Ip -|- 




Svrm^ 


/Vf2 

log SUSY 


+ 




mr 


/Vf2 


1 - 


xi 


12M|ugY/ j 


(4.1) 


where Xt = At — p/ tan/3 and Mg^gy = mq^rriq^, stays approximatively valid in the NMFV 
regime, so that peaks defined by \Xt\ ~ \/6Msusy are expected (see, e.g. Ref. [74] and 
references therein). 

Moving on with the slepton mass parameter M^, we observe a peak centred at around 
600 GeV after imposing all experimental constraints. This is mainly inferred by the anoma¬ 
lous magnetic moment of the muon requirement that strongly depends on the slepton sector 
properties. Turning to the Higgs sector (second line of Figure 1), the prior distribution of 
the /i-parameter shows a preference for low values while its posterior distribution slightly 
peaks around p ~ 200 GeV due to the Bg —)■ pp, Ao^ and AMb^ constraints which all de¬ 
pend on the chargino and neutralino sector. Next, the tan/3 parameter tends towards lower 
values both in its prior and posterior distributions, the favourite values being pushed to 
satisfy 12 < tan /3 < 18. Finally, the posterior distribution of the mass of the pseudoscalar 
Higgs boson rriA is shifted towards higher values with respect to its prior distribution. This 
results from the interplay of most considered observables for which low values of niA would 
induce too large Higgs contributions. 

The last two lines of Figure 1 concern the soft squark mass parameters. Low values 
are preferred for the first and second generation squark masses ^^ and ^, a 
feature that is mostly caused by the Higgs boson. This behaviour can be understood from 
the limiting case in which M% ~ M? ~ M? = m?. The one-loop corrections to nih 

VI,2 Ul,2 Ul,2 

that are proportional to 6^^ are here approximately given by [26] 


Ami = 




87r2(u^ 


+ vl) 


(^m)23 

(Yt^ 

(^«)23A 

ifi^ 

[ 2 

12m2 J 


(4.2) 


while the corresponding contributions of down-type squarks are obtained by replacing 
by Trf, Yt by 1), and by exchanging Vu and v^- In our parameterisation, 


(^«)23 - 


(4.3) 


so that for non-zero the Higgs boson becomes tachyonic if is too large. Similarly, 
the requirement of a physical solution for the electroweak vacuum also favours lower values 
for Mq^ ^. The distributions of the third-generation mass parameters AIq^ and prefer 
in contrast larger values due to AMb^ and the mass of the Higgs boson constraints. Finally, 
both the prior and posterior distributions of the right-handed down-type squark mass 
prefer lower values and are in this case very similar. 
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Figure 2. Same as Figure 1 in the case of the flavour-violating input parameters of our NMFV 
MSSM description. 

4.2 Flavour-violating parameters 

We now turn to the analysis of the constraints that are imposed on the seven non-minimally 
flavour-violating parameters 5^^ that are at the centre of interest of the present analysis. 
The corresponding prior and posterior distributions are displayed in Figure 2, and we detail 
the impact of the most important observables on Figure 3, Figure 4 and Figure 5. 

The theoretical constraints on any additional stop-scharm mixing in the left-left sector 
are relatively mild such that an almost flat behaviour is observed (see Figure 2). The 
5ll parameter is then mainly constrained by the i3-meson oscillation parameter AMb^ 
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Figure 3. Most relevant observables constraining the Sll (upper panel) and (lower panel) 
parameters. 


(which favours smaller absolute values of 5ll) and the branching ratio for the Bs — >• 
decay (which causes a slight preference to positive values), as shown in Figure 3. Values 
ranging up to |(5 ll| = 0.8 can nevertheless be reached, but this simultaneously requires large 
values for other 6 quantities so that cancellations between the different contributions to 
the considered observables occur (see Section 4.3). In a similar way, the prior distributions 
of the parameters <5^^ and <5^^ show a mild preference for low absolute values. The 
posterior distribution of the <5^^ parameter does not differ much from its prior distribution 
so that 5^^ is not sensitive to the experimental constraints under study. In contrast, 
the i?-meson oscillation parameter AMb^ restricts the posterior distribution of to be 
narrower while the Bg —)• /r/r branching ratio implies a preference to negative values (see 
Figure 3). However, the full explored range of — 0.8 < < 0.8 stays accessible in the 

context of both right-right mixing parameters. 

The flavour-violating left-right and right-left elements of the up-type squark mass 
matrix and 5'^b) turn out to be mainly constrained by the necessity to incorporate 
a Higgs boson with a mass of about 125 GeV, as can be seen in Figure 4. The posterior 
distribution of <5^^ exhibits two peaks at ~ 0.5 and is restricted to —0.15 < < 

0.15. Theoretically, this behaviour is expected from Eq. (4.3). The parameter however 
receives extra constraints stemming from the BR(i?s —)• observable (see Figure 4) so 


- 13 - 









Figure 4. Most relevant observables constraining the (left panel) and (centre and right 
panel) parameters. 


that the posterior distribution peaks around zero and has a maximal value of \Si^i\ ~ 0.2. 
We recall that the two parameters and (5^^ are independent and induce different 
mixing patterns. More precisely, describes a CL-tji mixing, while <5^^. corresponds to 
mixing between the cr and eigenstates. The impact of the constraints and the resulting 
distributions are therefore different and directly related to the structure of the chargino- 
squark-quark and neutralino squark-quark interactions. 

In the down-type squark sector, the prior distributions of the <5^^ and 6'^^ mixing 
parameters show a clear peak for values close to zero. Large values are often discarded as 
they imply large off-diagonal terms in the mass matrix so that the resulting spectrum 
likely contains tachyons. Both parameters are hardly constrained by any of the observables 
under consideration and we only observe minor effects. The posterior distribution of 
slightly prefers negative values, and the posterior distribution of is slightly narrower, 
when both distributions are compared to their respective prior. This mostly results from 
an interplay of all observables, although but for the 5^^ case, the B-meson oscillation 
observable AMr^ and the Higgs boson mass requirement play a non-negligible role (see 
Figure 5). 

We now illustrate the global distribution of all NMFV parameters. To this end, we 
introduce the quantities 


1^1 = [{SllY + + (4h)' + {SIr}" + + {5iRf + {SiRY 

login^l =log 6LL5lR5iR5lR5lR5iR5iR 


11/2 


(4.4) 


The former, |(5|, corresponds to the norm of a vector whose components are the seven NMFV 
parameters. Its value gives a measure of how far a given benchmark is situated from the 
minimally flavour-violating setup where |(5| =0. The maximum value that can be reached 
in our scan is |(I| ~ 1.56. The second quantity, loglH^I, corresponds to the logarithm of 
the absolute value of the product of the seven NMFV parameters. The case where all 
NMFV parameters are maximum corresponds to log jH^I ~ —3.5. In Figure 6, we show the 
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Figure 5. Most relevant observables constraining the parameter. 




Figure 6. The one-dimensional prior (yellow histogram) and posterior (violet curve) distributions 
of the quantities |i5| and loglll^l defined in Eq. (4.4). 


prior and posterior distributions of these two quantities. All scanned points feature |(5| > 0 
so that at least one of the NMFV parameters is sizeable and non-vanishing. The second 
quantity is in general large and negative so that at least one of the NMFV parameters has 
to be small. However, since the distribution shows a peak around log jH^j ~ —7 it is clear 
that a large fraction of the scanned points exhibit seven non-vanishing (with some sizeable) 
NMFV parameters. 

4.3 Correlations within the flavour-violating parameters 

Having discussed the distribution of single parameters, it is interesting to investigate pos¬ 
sible correlations between different NMFV quantities. A correlation indicator between two 
parameters x and y can be computed as 

, _ - x){yi - y) 

\/Er=i(^.-^)VEr=i(2/.-y)^ ’ 
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Figure 7. Two-dimensional distributions of the mostly correlated pairs of NMFV parameters after 
including all constraints. 


where the sum runs over all the points (xj, y*) of the sample and x and y are the mean values 
of the two parameters. The correlation factor is vanishing when there is no correlation, 
while r = ±1 indicates a linear correlation with the exception of the case in which the 
sampled points lie on a line parallel to one of the x and y axes. As the study of the 
correlations only makes sense when the parameter ranges cover the entire distribution 
spreads, we restrict our analysis to the NMFV parameters. The correlation indicators 
have been computed for any pair out of the seven NMFV parameters and the results are 
shown in Table 4 when one only accounts for the theoretical prior (second column) and 
after imposing the full set of constraints (last column). The correlations are found not 
particularly pronounced with all r-values being close to zero. 

We illustrate the correlations between different NMFV parameters on Figure 7. We 
however only focus on cases where the correlation indicator is above |r| > 0.25, namely 
on the This shows that scenarios in which 

several NMFV parameters are non-zero (and even significantly large) simultaneously are 
still allowed by current low-energy flavour and Higgs data. 

4.4 Squark masses and flavour decomposition 

We discuss in this section the distributions of the masses of the squarks, their flavour 
decomposition and the mass differences between states relevant for the LHC phenomenology 
of NMFV MSSM models. Figure 8 shows the prior and posterior distributions for the up- 
type squark masses. The shapes of the distributions for the two lightest states ui and U 2 
are very similar and they both peak at about 800-1000 GeV. The two lightest up-type 
states ui and U 2 , that are mostly of the first and second generation (see Figure 9), are in 
general relatively close in mass. This is due to the choice of common mass parameters for 
the first and second generation squarks. The heavier u-^, and states exhibit more 
spread distributions, the masses ranging from 1 to 3.5 TeV. Finally, the heaviest state 
uq is barely reachable at the LHC, with a mass lying in general above 2 TeV. Although 
the considered experimental constraints affect all NMFV supersymmetric parameters, the 
associated effects on the mass eigenvalues is at the end only mild, the mass distributions 
being only slightly shifted towards higher values. 
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Parameters 

th 

th -|- exp 

^Rl) 
i^RL^^R l) 
i^RR^^R l) 

-0.003 

0.007 

-0.000 

-0.002 

0.009 

0.270 

0.267 

-0.254 

0.185 

-0.158 

i^RR^ ^Rl) 
i^LL, 

i^RR’ ^Lr) 
i^LR^ ^Rl) 

i^LR^^L r) 
i^RR’ ^Lr) 
i^RR^ ^Rl) 

0.003 

0.002 

-0.021 

-0.001 

-0.002 

-0.024 

-0.006 

0.003 

-0.000 

-0.037 

-0.031 

-0.028 

0.027 

0.023 

0.013 

-0.012 

0.010 

-0.010 

i^RR^ ^Rl) 

i^RR’ ^Rr) 
i^RR^ ^Lr) 
i^RR^ ^Lr) 
i^LR^^R l) 

-0.002 

0.002 

0.001 

0.000 

-0.001 

0.000 

-0.008 

-0.007 

-0.006 

-0.003 

0.002 

0.000 


Table 4. The correlation coefficient r defined in Eq. (4.5) for all pairs of NMFV parameters. The 
parameter pairs are ordered by their correlation indicators when taking into account all imposed 
constraints (‘th+exp’). We also display the indicator values when only the theoretical prior is 
imposed (‘th’). 


From a phenomenological point of view, it is interesting to examine the flavour (in 
particular the stop) content of the six up-type squarks. The posterior distribution of 
the stop content of the three lightest up-type squarks is depicted in Figure 9 and shown 
in correlation with the respective squark mass. The lighter states ui, U 2 -, us (and also 
U 4 ) are mainly not stop-like, i.e. they have a significant up or charm component. Most 
scanned scenarios indeed exhibit a charm-dominated lightest ui squark, while U 2 is mostly 
dominated by its up component. This contrasts with usual flavour-conserving MSSM 
setups where the lightest squark state is typically a stop. This feature can be traced to 
the first and second generation soft masses that are driven to lower values as explained in 
Section 4.1, whilst the third generation squark masses are pushed towards higher values 
by the flavour constraints. Furthermore, even in the presence of large trilinear terms, the 
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Figure 8. One-dimensional prior (yellow histogram) and posterior (violet curve) distributions of 
the masses of the six up-type squarks. 




Figure 9. Resulting correlations between the stop flavour content and the masses of the three 
lightest up-type squarks after imposing all experimental constraints mentioned in Table 3. Red 
colour indicates the highest and dark purple the lowest likelihood. 


lightest states are still found to be up-like or charm-like. 

Similar conclusions hold for the sector of the down-type squarks. We show their masses 
in Figure 10 and selected flavour decompositions in Figure 11. The three lighter states ex¬ 
hibit comparable distributions, peaking as for the up-type squarks at about 800-1000 GeV. 
The mass distributions of the and states feature distributions with a larger spread, 
and the one of heaviest dg squark is peaking at about 3 TeV, although masses of about 
1 TeV are predicted for a small subset of scenarios. Flavour mixing in the down-type squark 
sector is generally less pronounced than for the up-type squarks, as illustrated on Figure 11 
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Figure 10. Same as Figure 8 for the down-type squark sector. 




Figure 11. Resulting correlations between the sbottom flavour content and the masses of the three 
lightest down-type squarks after imposing all experimental constraints mentioned in Table 3. Red 
colour indicates the highest and dark purple the lowest likelihood. 


where we depict the correlations between the sbottom content and the masses of the lighter 
down-type squarks. A majority of scenarios include light down-like and strange-like squark 
states and there is only a small number of parameter points where di and d 2 contain a 
sizeable sbottom content. 

In Figure 12, we show the correlations between the masses of the lightest squark states 
and the one of the lightest neutralino. For most (~ 95%) viable points, the mass difference 
is well above 50 GeV, which is a favorable condition for collider searches as the spectrum is 
not compressed. A considerable number (~ 40%) of parameter points features ui masses 
of about 500 - 1000 GeV together with neutralino masses of the order of 150 - 400 GeV. 
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Figure 12. Correlations between the lightest neutralino mass and the lightest up-and down-type 
squark masses. Red colour indicates the highest likelihood. 


Such mass configurations are likely to be ruled out by Run I LHC data. This is accounted 
for in the next section, where we include collider constraints on the NMFV MSSM setup 
and define benchmark scenarios suitable for searches at the LHC Run H. 


5 Benchmark scenarios 

In this section, we identify an ensemble of benchmark scenarios capturing typical features of 
the parameter space regions favoured by the constraints previously investigated. We have 
ordered all acceptable parameter setups according to their likelihood and selected four 
scenarios among the best ones. Our selection is aimed to cover different phenomenological 
properties of the NMFV MSSM and to be relevant for future LHC searches. The input 
parameters corresponding to the benchmark scenarios of our choice are indicated in Table 5. 
In Figure 13, we present the mass spectra of the four selected scenarios and depict the 
flavour content of the different squark eigenstates. We hnally show in Table 6 the branching 
ratios related to the dominant decay modes of the squarks ligther than about 1 TeV. 
Additionally, we have verified that the electroweak vacuum is stable for all selected points 
by using the programme Vevacious [77]. We now briefly outline the main characteristics 
of the four proposed benchmark scenarios^. 

Scenario I This benchmark point presents one up-type and two down-type squarks with 
masses below 1 TeV. The lightest up-type squark is mostly stop-like, although it contains 
a small scharm component, and has a mass of 831 GeV which implies a sizeable production 
cross section at the LHC. In the sector of the down-type squarks, the two lightest states 
are almost purely sbottom-like with masses of 763 GeV and 854 GeV respectively. Since 
only the heaviest neutralino and chargino states are heavier than the three lightest squark 
states, various decay channels are open so that the real challenge for future LHC analyses 
would be to become sensitive to flavour-violating branching ratios of a few percents. In 

^The benchmark scenarios can be provided, under the form of Supersymmetry Les Houches Accord 
compliant files [78], by the authors upon request. 
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Parameter 

I 

II 

III 

IV 

as{mz) 

1.187- 10“^ 

1.194-10-1 

1.176 - 10-1 

1.176 - 10-1 

pole 

in'- 

176.00 GeV 

175.53 GeV 

173.53 GeV 

174.02 GeV 

mb{mb) 

4.10 GeV 

4.24 GeV 

4.28 GeV 

4.10 GeV 


1192.7 GeV 

2288.2 GeV 

637.7 GeV 

753.2 GeV 


883.7 GeV 

425.3 GeV 

3483.0 GeV 

2662.7 GeV 


2412.6 GeV 

1757.7 GeV 

934.0 GeV 

984.7 GeV 

^Pa 

2344.3 GeV 

2753.8 GeV 

2862.2 GeV 

2010.6 GeV 


2295.1 GeV 

551.6 GeV 

1331.1 GeV 

882.7 GeV 

^Pa 

843.8 GeV 

713.5 GeV 

901.8 GeV 

670.5 GeV 


-2424.1 GeV 

1807.3 GeV 

1586.3 GeV 

-2833.4 GeV 

tan j3 

17.4 

21.1 

29.2 

34.0 

P 

615.7 GeV 

772.8 GeV 

508.1 GeV 

442.7 GeV 

rriA 

1334.5 GeV 

1300.3 GeV 

1294.8 GeV 

1431.0 GeV 

Ml 

474.5 GeV 

315.3 GeV 

525.2 GeV 

390.0 GeV 

Me 

2466.5 GeV 

1552.5 GeV 

3396.7 GeV 

2813.4 GeV 

Sll 

1.4- 10"^ 

-4.6 - 10-2 

3.7-10-1 

5.9 - 10-1 

^RR 

1.7-10-1 

2.2 - 10-1 

7.3 - 10-1 

6.0 - 10-1 

Xd 

^RR 

1.4-10-1 

-1.4-10-1 

-2.9 - 10-1 

-7.5 - 10-1 

^Ir 

9.2 - 10-2 

3.5 - 10-2 

1.7-10-1 

1.0 - 10-1 

^LR 

-3.9 - 10-2 

-3.6 - 10-3 

-6.1 - 10-3 

4.6 - 10-3 


-9.7-10-2 

1.4-10-2 

-9.9 - 10-2 

-7.2 - 10-2 

^RL 

-7.6 - 10-1 

-1.4-10-2 

-1.2 - 10-3 

1.6 - 10-3 


Table 5. Definition of four benchmark points suitable for phenomenological studies of the NMFV 
MSSM. 


this scenario, the electroweak vacuum is long-lived and has a lifetime larger than the age 
of the Universe. 

Scenario II In this scenario, only the lightest of the up-type squarks is expected to lie 
within the reach of LHC, with a mass of 526 GeV. It is almost a pure stop state with a 
small charm component. Since the only lighter superpartner is the lightest neutralino, it 
will preferably decay into a system. There are four down-type squarks lying below 
1 TeV, their masses being 519 GeV, 555 GeV, 566 GeV and 747 GeV. These four states 
are admixtures of all three flavours and their dominant decay modes include in particular 
final states containing the next-to-lightest neutralino or the lightest chargino. Contrary 
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Figure 13. Mass spectra of the benchmark scenarios defined in Table 5. The colour code that has 
been employed for depicting the squark eigenstates indicates their flavour content: the {red, green, 
blue} colour corresponds to the {first, second, third} generation flavours. 

to the scenario I, the branching ratios related to flavour-violating decays can reach up 
to 16 percents, which make them possibly testable at the LHC. Moreover, the fourth 
down-type squark has sizeable branching ratios for decays into the lightest up-type squark 
and a VL-boson as well as into the lightest down-type squark and either a Z-boson or a 
Higgs boson. Although many squarks are very light, this scenario evades all LHC Run I 
constraints thanks to a heavy lightest neutralino of 315 GeV. Moreover, the vacuum has 
been found to be stable. 

Scenario III This benchmark point features numerous squark mass eigenstates in the 
reach of the LHC. It indeed exhibits three up-type squarks with masses of 836 GeV, 
882 GeV and 928 GeV, and three down-type squarks with masses of 880 GeV, 931 GeV 
and 1050 GeV. The stop-like states are here the heaviest ones, and the up-type squark 
states reachable at LHC only contain up and charm flavours. Similarly, the heavier and 
dg down-type squarks are the only ones containing a sbottom component. This feature 
is the direct consequence of the lower values that are favoured for the ^^ and 
^ parameters, when compared to the values that are favoured by the third-generation 
soft parameters Mq^, and Equivalently, this can be seen as an implication 

of allowing for flavour-violating entries in the squark mass matrices (see Section 4.1). In 
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Decay 

I 

II 

III 

IV 

ui txi 

0.14 

0.99 

0.06 

0.09 

ui cxi 


0.01 

0.24 

0.07 

ui tx^ 

0.24 


0.11 


ui txl 

0.43 


0.06 

0.26 

ui tx^ 




0.07 

ui cxs 



0.27 

0.10 

-t bxt 

0.19 


0.24 

0.22 

ui -t 



0.02 

0.17 

ui —>• W^di 




0.02 

U2 -t tXi 



0.07 

0.04 

U2 -t CXi 



0.22 

0.32 

U2 tx^ 



0.08 

0.09 

U2 txi 



0.09 

0.11 

to 

coo 



0.35 

0.06 

U2 -t bxt 



0.11 

0.21 

U2 -t 



0.09 


U2 — t W^d2 




0.06 

U2 — )• Z'^Ui 




0.05 

U2 —t h^Ul 




0.02 

U3 -t CXi 



0.03 

0.09 

U3 -t UXi 



0.03 

0.03 

U3 -t txi 




0.02 

U3 txi 




0.02 

U3 -t cxi 



0.05 

0.02 

U3 -t uxl 



0.45 

0.09 

U3 -t CX4 




0.01 

U3 -t uxl 




0.15 

U3 -t bxt 



0.01 

0.05 

U3 -t dxt 



0.41 

0.01 

U3 -t dxt 




0.33 

U3 — >■ W^d2 




0.01 

U3 —)■ Z^UI 




0.02 


Decay 

I 

II 

III 

IV 

di bxl 

0.73 

0.84 


0.35 

di -t sXi 


0.16 

0.40 

0.09 

di bxl 

0.15 


0.01 

0.30 

di SX2 



0.02 


di bxl 

0.12 



0.25 

di sxl 



0.05 


di -t txr 



0.42 


di -t cXi 



0.09 


d2 bxl 

0.04 

0.05 


0.02 

d2 -t sxi 




0.24 

d2 -t dxi 



0.70 


d2 bxl 




0.04 

d2 -t SX2 

0.04 

0.95 



d2 dxl 



0.03 


d2 dxl 



0.08 


d2 bxl 




0.04 

d2 sxl 

0.04 



0.02 

d2 -t txi 

0.87 



0.51 

d2 -t cXi 




0.09 

d2 -t uxi 



0.15 


d2 — t IV“i6i 



0.03 


d2 Z^di 




0.02 

d2 h^di 




0.02 

d3 -t dXi 


1.00 


0.16 

d3 -t dxl 




0.01 

d3 -t dxl 




0.25 

d3 -t uxi 




0.07 

ds -t CX2 




0.02 

d3 -t UX2 




0.48 

di bxi 


0.42 



di -t sXi 


0.03 



di dxi 




0.83 

di -t dx3 




0.17 

di —>■ W~ui 


0.27 



di Z^di 


0.13 



— )* hP d\ 


0.13 



d4. —)* d2 


0.01 




Table 6. Branching ratios associated with the dominant decay modes of the up-type (left) and 
down-type (right) squarks lighter than about 1 TeV for the benchmark points defined in Table 5. 
Branching ratios below 1% are not indicated. 
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addition, all gauginos except the heaviest neutralino and chargino feature lower masses, 
so that a variety of decay channels are open. Finally, we found a direction in which, for 
extremely large field excitations, the electroweak vacuum is unbounded from below. This 
situation is similar to the case of the Standard Model [79, 80], and higher order corrections 
to the scalar potential would be needed to make a conclusive statement. 

Scenario IV Our last scenario features numerous squark states as well as a complete 
electroweakino spectrum below 1 TeV. More precisely, the lighter up-type squarks have 
masses of 751 GeV, 902 GeV and 923 GeV, while the lighter down-type squark masses are 
of 582 GeV, 775 GeV and 924 GeV. In addition, three other states are not too far above 
1 TeV with masses of 1119 GeV (ua), 1029 GeV (^ 4 ) and 1167 GeV (ds). The two lightest 
up-type squarks consist in this case of a mixture of all three flavours (with a dominant 
charm content), which leads to interesting decay patterns, as shown in Table 6 . Similarly 
to scenario I, this scenario exhibits a long-lived electroweak vacuum with a lifetime larger 
than the age of the Universe. 

6 Conclusion 

We have studied non-minimal flavour-violation in the MSSM by allowing for flavour mixing 
between the second and third generation squarks. We have used a Markov Chain Monte 
Carlo scanning technique to explore the underlying parameter space and imposed a set of 
experimental constraints arising from U-meson and kaon physics. We have additionally 
enforced the model to accommodate a light Higgs boson with a mass of 125 GeV. 

First and second generation soft squark masses are theoretically restricted to low values 
in order to avoid tachyons in the Higgs sector. As a consequence, the lighter squarks 
are often not the stop and sbottom ones, which contrasts with scenarios of the usual 
minimally flavour-violating MSSM. Requiring a theoretically consistent Higgs sector and a 
light Higgs boson of about 125 GeV similarly restrict the left-right and right-left flavour- 
violating squark mixing parameters to be small. In contrast, the 6ll and 

NMFV parameters are mainly constrained by neutral H-meson oscillations, and the rare 
Bg —)• /i/i decay mainly influences All other NMFV parameters are left unconstrained 
by the considered experimental observations. 

In view of the recently started second LHC run, we have used our MCMC scan results 
to propose four benchmark scenarios allowed by current data that exhibit distinct features 
and that are suitable for future analyses of NMFV effects in the MSSM. In most proposed 
scenarios, several squarks have masses close to 1 TeV so that they should be reachable 
within the next few years. 
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